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The impact of magnetic field on the cluster M-T relation 
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Abstract. We discuss the impact of magnetic field on the mass - temperature relation for groups and clusters 
of galaxies based on the derivation of the general Magnetic Virial Theorem. The presence of a magnetic field B 
yields a decrease of the virial temperature T for a fixed mass M: such a decrease in T is stronger for low-mass 
systems than for high-mass systems. We outline several implications of the presence of B-field and of its mass 
scaling for the structure and evolution of groups and clusters. 
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1. Introduction 

Magnetic fields fill intracluster and interstellar space, af- 
fect the evolution of galaxies, contribute significantly to 
the total pressure of interstellar gas, are essential for the 
onset of star formation, and control the diffusion, the con- 
finement and the evolution of cosmic rays in the interstel- 
lar and intracluster medium (ICM). In clusters of galaxies, 
magnetic fields may play also a critical role in regulating 
heat conduction (e.g., Chandran et al. 1998, Narayan & 
Medvedev 2001), and may also govern and trace cluster 
formation and evolution. 

We know that magnetic fields exist in clusters of galax- 
ies for several reasons. First, in many galaxy clusters we 
observe the synchrotron radio-halo emission produced by 
relativistic electrons spiraling along magnetic field lines. 
Second, the Faraday rotation of linearly polarized radio 
emission traversing the ICM proves directly and indepen- 
dently the existence of intracluster magnetic fields (see, 
e.g., Carilli & Taylor 2002, Govoni & Feretti 2004 for re- 
cent reviews). The Rotation Measure (RM) data through- 
out the inner (~ 0.5 Mpc) cluster region support magnetic 
field strengths of the order of several to tens of /j,G (see 
Carilli & Taylor 2002, Govoni & Feretti 2004). The high 
local values of B observed in the central, cool region of 
clusters are likely related, however, to quite special con- 
ditions (such as turbulent amplification of the local B- 
field driven by radio bubbles or AGN jets, see e.g. Ensslin 
& Vogt 2006) and thus are probably not representative 
of the overall system (see, e.g. Carilli & Taylor 2002). 
Other estimates of the magnetic field strength on the clus- 
ter wide scale come from the combination of synchrotron 
radio and inverse Compton detections in the hard X-rays 
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(e.g., Colafrancesco, Marchegiani & Perola 2005), from the 
study of cold fronts and from numerical simulations (see, 
e.g., Govoni & Feretti 2004). This evidence provides indi- 
cation on the wide-scale B-field which is at the level of a 
few tens up to several fiG (and in some cases up to ~ 10 
/iG, as in Coma) with the larger values being attained by 
the most massive systems. 

Numerical simulations (e.g., Dolag et al. 2001a) have 
shown that the wide-scale magnetic fields in massive clus- 
ters produce variations of the cluster mass at the level 
of ~ 5 — 10% of their unmagnetized value. Such mass 
variations induce a comparable variation on the IC gas 
temperature T for virialized systems. Such variations are 
not expected to produce strong variations in the relative 
M — T relation for massive clusters. 

The M — T relation predicted in a pure CDM model 
for B = follows the self-similar scaling M oc T v with 
r\ = 3/2 (see, e.g., Colafrancesco et al. 1997, Arnaud 
2005). A Chandra study (Allen et al. 2001) of five hot 
clusters (with ksT g > 5.5 keV) derived a M — T relation 
slope of -q — 1.51 ± 0.27, consistent with the self-similar 
model. However, due to the relatively small Chandra field 
of view, the M — T relation was established at -R25001 
i.e., about 0.3i?20o (here Rs and M§ are the radius and 
mass at which the density contrast of the system is S). 
More recently, the M — T relation was established down 
to lower density contrasts (6 = 200) from a sample of 
ten nearby relaxed galaxy clusters covering a wider tem- 
perature range, fceT g w 2 — 9 keV (Arnaud et al. 2005). 
The masses were derived from mass profiles measured with 
XMM-Newton at least down to -R1000 and extrapolated be- 
yond that radius using the NFW (Navarro, Frenk & White 
1997) model. The M2500 — T for hot clusters is consistent 
with the Chandra results. The slope of the M — T relation 
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is the same at all 5 values, reflecting the self-similarity of 
the mass profiles. At S = 500 the slope of the relation 
for the sub-sample of hot clusters (ksTg > 3.5 keV) is 
rj = 1.49 ± 0.15 consistent with the standard CDM self- 
similar expectation. The relation, however, steepens when 
the whole sample of clusters is considered, providing a 
slope rj = 1.71 ± 0.09. The normalisation of the M — T re- 
lation differs, at all density contrasts from the prediction 
of pure gravitation based models by ~ 30% (see Arnaud 
2005 for a discussion). 

In this Letter we will explore the effect of wide-scale mag- 
netic fields on the M — T relation over a large range 
of masses and temperatures by using the predictions of 
the magnetic virial theorem. We will discuss its implica- 
tions for the evolution and the scaling relations of magne- 
tized clusters. The relevant physical quantities are calcu- 
lated using Hq = 71 km s _1 Mpc -1 and a flat, vacuum- 
dominated CDM (S! m = 0.3, £1a = 0.7) cosmological 
model. 



where 1(c) = J^(p g (r = 0)/p g (z = 0)) 2a x 2 y 2 g a (*, B = 
0)dx. Here c = r V i r /r s (we assume a NFW Dark Matter 
density profile with scale radius r s ) and y g (x,B = 0) = 
p g (x) / p g (x — 0) is the gas density profile normalized to 
the central gas density (i.e. the solution of the hydrostatic 
equilibrium equation in the absence of magnetic field, see 
Colafrancesco & Giordano 2006a for details). The radial 
profile of the magnetic field has been assumed as B(r) = 
B if [p g (r,B)/\Q A p(z = 0)]" with a = 0.9 (see, e.g., Dolag 
ct al. 2001b). 

For the case P ex t — and B = 0, the quantity M$ = 
and the well-known relation 
£pm p W 



k B T g (B = 0) 



(5) 



re-obtains (here p — 0.63 is the mean molecular weight, 
corresponding to a hydrogen mass fraction of 0.69, m p is 
the proton mass and fee is the Boltzmann constant). 
For B > 0, the quantity Mq > and the gas temperature 
at fixed M v i r , as obtained from eq.©, is 



2. The magnetic virial theorem for galaxy clusters kT g — kT g (B — 0) 1 
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Under the assumption of a ICM in hydrostatic equilib- 
rium with the potential well of a spherically-symmetric, 
isolated, virialized and magnetized cluster, the general re- 
lation between the ICM temperature T and the cluster 
virial mass M is obtained by applying the magnetic virial 
theorem (MVT): 



1 d 2 L, 



2 dt 2 
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ik 



F ik d 3 x + W ife , 



(1) 



where 1^ is the inertia momentum tensor, Kik is the ki- 
netic energy tensor, U is the thermal energy of the intra- 
cluster gas, is the Maxwell tensor associated to the 
magnetic field and Wiu is the potential energy tensor. The 
full derivation of the MVT is reported in the Appendix. 
For a static and isothermal galaxy cluster the trace of eq. 
(fT]) yields the condition 



2K + 2U + U B + W = . 



(2) 



where Ub is the magnetic energy of the system, U is the 
kinetic energy of the gas, K is the dark-matter particle 
kinetic energy and W is the potential energy of the system 
(see Appendix for details). The previous eqs. © and (J2j) 
hold specifically in the absence of an external medium. 
For the general case of a cluster which is immersed in a 
Inter Galactic Medium (IGM) or external medium which 
exerts an external pressure P ex t > eq. ((2]) yields the formula 
for the temperature of the gas in virial equilibrium 
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and is lower (for P ext = 0) than that given by eq.© be- 
cause the additional magnetic field energy term U b adds 
to the MVT. The presence of an external pressure P ex t 
tends to compensate the decrease of T g induced by the 
magnetic field. For values of the temperature and density 
of the IGM (as estimated by the WHIM structure around 
large-scale overdensities, see, e.g., Fang & Bryan 2001), 
P ext ~ 1.7 • 10" 3 eV cm- 3 (n IGM /10- 5 cm- 3 )(T IGM /2 ■ 
10 6 K). However, in the outer regions of massive clusters 
(at r < r V i r ) the external gas pressure can reach values 
P ext - 0.2 eV cm- 3 (n/10- 4 cm- 3 )(T a /1.7 • lO 7 ^) [here 
we considered the mean projected temperature profile for 
the cluster sample studied by Piffaretti et al. (2005, see 
their Fig. 4) and a typical cluster with Tx = 10 keV]. 
In such a case, the value of P ex t is a significant fraction 
- 4% of the central ICM pressure and - 50% of the ICM 
pressure at the virial radius for a typical cluster. Thus, 
it cannot be neglected in the T g estimate from eq.©. A 
value P ext ~ 0.2 eV cm -3 , as estimated at the outskirts 
(r ^ r v i r ) of rich clusters, can be considered as an upper 
bound to P ex ti since an exact determination of the total 
cluster mass (which is subject to various systematic un- 
certainties, see, e.g., Rasia et al. 2006) certainly requires 
to go beyond r vir . We thus consider in the following this 
value of P ex t as a reference upper bound to be used in 
our temperature estimate from eq.© in the presence of 
a B-field. Lower values of P ex t, down to its value in the 
WHIM, have progressively minor importance. 
For reasonable values of B* ^ a few pG, the quantity 
M 2 > M 2 mr ■ (P ext /P mr ) (with P mr = (ig^)- 1 ) in 
eq.©, and the main effect is a reduction of the cluster 
temperature which is more pronounced for less massive 
systems, where M^ becomes comparable to M vir . The ef- 
fect of the magnetic field and of the external pressure are 
larger for low-M clusters (see Fig©. 
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3. The magnetized M-T relation 

The T — M relation for magnetized clusters is shown in 
Fig[TJ We normalize the T spec t r . — M200 relation for the 



10 



0.1 







1 1 




1 




























































/ 1 








/ 1 




/ 




i 1 




-5 / 




1 

i 




" / 


i 


1 




- / 


i 

i 

i 


I 




■h 




i 








i 
i 




7 










j20 








j 


30 

1 1 





I — — , — — I — I — — I 1 — 

0.1 1 10 



M aoo /M B 

Fig. 1. We show the T apec t r — M200 relation at z = for 
clusters which contain a magnetic field _B„ in the illustrative 
range — 30 [iG (as labelled) in the case of P ext = 0.2 eV 
cm " 3 . Here M 8 ~ 2-10 14 M Q fcf 1 1 . Data are taken from Arnaud 
(2005). 



case B = to the observed data derived by Arnaud et al. 
(2005) by assuming T spectr = T g with £ w 1.5 in eq.©, 
as can be expected from the continuous shock-heating of 
the 1C gas within the virial radius after the formation of 
the original structure (see, e.g., Makino et al. 1998, Fujita 
et al. 2003, Ryu et al. 2003). Such a prescription for £ 
has been used by the previous authors for unmagnetized 
clusters in the absence of external pressure. A systematic 
effect which goes towards the direction of increasing the 
cluster temperatures is the presence of a minimal external 
pressure in eq.®. A value P ext ~ (0.1 — 0.2)P„ r (like that 
found in the IGM around clusters) could easily accommo- 
date for an overall value of £ « 1.5 x (P ex t/P V ir) > 1.5 
and reasonably in the range 1.65 — 1.8, for the previous 
values of P ex t- Given the large theoretical uncertainty on 
the non-gravitational heating efficiency, we adopt an over- 
all value £ « 1.8 to normalize our prediction to the data 
point at M200 ~ 3- Ms in Fig.l, which is the point with the 
smaller intrinsic error. Values of £ in the plausible range 
1.5 — 1.8 marginally change, however, our predictions. 
The relation M200 — 0.77 M V i r is also found in our mass 
scale definition. 

Small variations of temperatures with respect to their un- 
magnetized values are found for massive clusters since the 
quantity <C M V i r in this mass range and the value 
of P ext has little or negligible effect (see Figj2]). This is 
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Fig. 2. Same as Fig[T]but for a P ext = with B, = (solid 
curves) and B* = 20 /iG (long-dashes curves) and P ex t = 0.2 
eV cm -3 for the same values of B„ = (short-dashes curves) 
and B* = 20 /iG (dotted curves) The evolution of the T apectr — 
M200 relation for magnetized clusters is shown at z = (black), 
z = 0.5 (red) and z = l(green). 



in agreement with the results of numerical simulations 
(Dolag et al. 2001a). However, when becomes com- 
parable to M v i r , the IC gas temperature becomes lower 
than its unmagnetized value and the T — M relation steep- 
ens in the range of less massive systems like groups and 
poor clusters. The temperature T„ formally tends to zero 
when — » M vir (l + P ex t/ P V ir) ■ However, this limit is 
unphysical since it corresponds to an unstable system in 
which the magnetic pressure overcomes the gravitational 
pull. Thus, any physical configuration of magnetized viri- 
alized structures must have < M„i r (l + P ex t/ Pvir) 1 ^ 2 ■ 
The effect of P ex t counterbalances the effect of the i?-field 
on the T — M relation, increases for low-M systems and 
decreases with increasing redshift (see Fig(2J) because P V i r 
increases with increasing redshift. 

4. Discussion and conclusions 

We have derived here, for the first time, a relation between 
the temperature of the IC gas from the general MVT in 
the presence of magnetic field and external pressure. The 
result of the MVT for clusters bring relevant modifications 
to the gas temperature for virialized and magnetized clus- 
ters. As a consequence, the observed T—M relation is 
steeper than the simple predictions of a ACDM structure 
formation scenario and its effective slope increases in the 
low-M region. However, since the masses of the observed 
clusters have been derived under the assumption of ab- 
sence of B field, the slope indicated by the data of the 
observed T — M relation could be not completely repre- 



4 



S. Colafrancesco & F. Giordano: Magnetic field and the cluster M-T relation 



sentative. In this context we also stress that the T — M 
relation might be affected by other systematic uncertain- 
ties in the mass derived by X-ray observations (e.g., Rasia 
et al. 2006) which would change the slope of the T — M 
relation especially in the low-M range. A robust analysis 
of the cluster mass estimate should require the use of a 
detailed hydrostatic equilibrium condition in combination 
with reliable temperature profiles. In both these aspects 
the effect of the B-ficld is relevant and should be taken into 
account. Furthermore, the predictions at low-T, where the 
effects of the B-field are stronger, are rendered uncertain 
by the absence of a clear definition of a spectroscopic tem- 
perature (e.g., Mazzotta et al. 2004). Thus, a complete 
analysis of the T — M relation relies on a very detailed 
understanding of the physical properties of the IC gas in 
the presence of B-field with the input of a precise total 
mass reconstruction and temperature determination. 

The results we derived here have a broad range of im- 
plications on cluster structure and evolution: flattening of 
the entropy - temperature relation and higher entropies 
in cluster cores are expected in the presence of magnetic 
fields. Further effects on the X-ray luminosity - tempera- 
ture relation are also expected as well as modifications of 
the thermal Sunyaev-Zel'dovich effect. Since these studies 
are far beyond the scope of this paper, we refer the in- 
terested reader to much more detailed analysis which are 
presented elsewhere (Colafrancesco & Giordano 2006a, b). 

To conclude, we notice that a full description of the 
structure and evolution of the population of groups and 
clusters of galaxies which considers also the role of mag- 
netic fields will definitely shed light on several, still unclear 
aspects of the interference between gravitational and non- 
gravitational mechanisms in the evolution of these sys- 
tems, and calls for a more refined physical description to 
use galaxy clusters as appropriate cosmological probes. 
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tract No.2004027755_003. 

Appendix A: The MVT for galaxy clusters 

Let us introduce the following quantities: 



ip(r) 



-G 



V \xj - Xj>\ 



d 3 x' 



••• - -H / f^f 

z Jy JV \ x 3 x 3'\ 



3 —d 3 xd 3 x' 



(A.l) 
(A.2) 



Ft 



8^ 



47T 



where ip is the gravitational potential, W is the gravita- 
tional energy, Fy is the Maxwell tensor associated to the 
magnetic field and p = p ( i m + Pg is the total density of the 
cluster (we neglect here the subdominant contribution of 
galaxies). Then, the equation of motion for the systems 
given by the Euler equation writes as 



„(»+(„. v)),, = -g-§ 



where Vi is the i-th component of the velocity and p is the 
total pressure given by pdm + Pg- Here, p ~ p g since we 
assume that DM is cold and collisionlcss, i.e. pd m ~ (we 
consider a fluid with no viscosity for which = p Sij). 
Multiplying by x k and integrating over the cluster volume 
we obtain: 

x k ^-{pv i )d 3 x 
v ot 

/3 f dp 

x k T^-{pviVj)d 3 x = - J x fe — d 3 x (A. 5) 

dip 3 f dF i:j 3 
x k p— — ax — / Xk ax. 



dxi 



dxj 



Using the continuity equation and the standard integral 
theorems, we convert the first member of this equation in 
the form 



d f d 

x k — (pv l )d 3 x + / x k - — (pviVj)d 3 x 
v at Jy oxj 

d 



(A.6) 



dt 



v 



pviVjd 3 x — 2K t k + (t> XkpviVjdSj 



where Kij = 1/2 J pd m ViVjd 3 x indicates the dark-matter 
kinetic energy tensor. The second member of eq.(A.5) 
writes as 



dp q 2 
Xh ~dx~ d X = 3 USik ~ t dSjVkP 



(A.7) 



OF 



J F lk d 3 x- j> dS-jF.jXk (A. 



- / x k —^d 3 x = I F lk d A x 

Jy OX j 

where U — 3/2 J p g d 3 x is the IC gas thermal energy and 
one can show that 



XkPj—d x = -Wik- 

y OXi I 



(A.9) 



Neglecting (in our case) the surface integrals (these phys- 
ical quantities are negligible when the integration surface 
is chosen far from the cluster center) one obtains: 



1 d 2 I lk 

2 dt 2 



= 2K lk + ~U8 ik + / F lk d 3 x + W lk , (A.10) 



where is defined as 

f 3 



(A.ll) 



Using the trace of eq. (|A.10[) . we obtain the equation 
ld 2 I 



(A.3) 2 dt 2 



2K + 2U + U B + W 



(A.12) 



where Ub = F, and F = J d 3 x. For a cluster in a 



static configuration (quite a good approximation for real 
systems) one has 



2 dt 2 ' 



(A.13) 



from which eq.([2) derives. 

In a more general derivation of the MVT (i.e., taking into 
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account also surface integrals in eqs.A7 and A8), the mag- 
netic energy writes as 



Ub - <t> XkF.jdSj = — , 



(A.14) 



where <f) = n(B*/ iiG)r% ir is the magnetic flux through the 
equatorial section of the system. The trace of the surface 
integral in ea. (|A.7[) writes as 



Pext <j> (r ■ dS) = P ext ^rl„ , 



(A.15) 



and is usually considered in the analysis of the standard 
Virial Theorem without the influence of a B-field (see, e.g., 
Carlberg et al. 1997). 
In the case of isothermal systems: 



2U = 3 / p g d 3 x = 3 c 2 s p g d 3 x ~ 3c 2 M s 
Jv Jv 



where 



c 2 _ k B T g 



(A.16) 



(A.17) 



The dark-matter particle kinetic energy writes as 

2K = {vjj [ Pdm d 3 x = {vj m )M dm , (A.18) 



where (w^) 1 ^ 2 is the dark-matter particle velocity disper- 
sion. From eqs. (A.12)-(A.15) we obtain: 



M 9 + (v 2 d JM dm = 



^ (l 



M- 



-P K 



Since l/VSiv^) 1 / 2 ~ c s , one obtains 



(A.19) 



^(Mg + M dm ) 



Air. 



•)■ 



from which eq.Q follows setting M g + M dv 



(A.20) 
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